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Abstract 

According to the Wulff construction the shape of the equilibrium crystal 
is determined by the surface tension considered as a function of the interface 
orientation. We present some (conjectured) approximate solutions and some 
rigorous results concerning this function, in the case of a lattice gas, and 
apply them to study the shape of the equilibrium crystal and, in particular, 
the shape of the facets of this crystal. 
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1 Introduction 

In a first approximation one can model the interatomic forces in a crystal 
by a lattice gas. In a typical two-phase equilibrium state there is, in these 
systems, a dense component, which can be identified as the crystal phase, 
and a dilute component, which can be identified as the vapor phase. The 
underlying lattice structure implies that the crystal phase is anisotropic, 
while this assumption, though unrealistic for the vapor phase, should be 
immaterial for the description of the crystal- vapor interface [T], [2]. As an 
illustrative example of such systems, the ferromagnetic Ising model will be 
considered. 

The Ising model is defined on the d- dimensional cubic lattice C = Z d } 
with configuration space Q = {—1, 1} C - The value a(i) = ±1 is the spin 
at the site i G C. The occupation numbers n(i) = (l/2)(cr(z) + 1) = or 
1 give the lattice gas version of this model. The energy of a configuration 
°a = £ A}, in a finite box A C £, under the boundary conditions 

a G Q, is 

H A (a A \a) = -J W 

where are pairs of nearest neighbor sites and a(i) = a(i) if i ^ A. The 
partition function, at the inverse temperature /3 = 1/kT, is given by 

Z°(K) = 5>xp ( - f3H A (a A | a)) (2) 

CTA 

The limit 

/(/?)= lim -J-ln^(A) (3) 

is independent of the boundary conditions and defines the free energy per 
unit volume. 

This model presents, at low temperatures T < T c , where T c is the criti- 
cal temperature, two distinct thermodynamic pure phases. This means two 
extremal translation invariant Gibbs states, which correspond to the limits, 
when A — > oo, of the finite volume Gibbs measures 

Z~ a {K)~ 1 e W {-(3H A {a A \ a)) (4) 

with boundary conditions a respectively equal to the ground configurations 
(+) and (— ) (respectively, a(i) = 1 and a(i) = —1, for all i G C). On the 
other side, if T > T c , the Gibbs state is unique. 
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Each configuration inside A can be described in a geometric way by spec- 
ifying the set of Peierls contours which indicate the boundaries between the 
regions of spin 1 and the regions of spin —1. Unit square surfaces are placed 
midway between the pairs of nearest-neighbor sites i and j, perpendicularly 
to these bonds, whenever a(i)a(j) = — 1. The connected components of this 
set are the Peierls contours. Under the boundary conditions (+) and (— ), 
the contours form a set of closed surfaces. They can be viewed as defects, 
or excitations, with respect to the ground states of the system (the constant 
configurations 1 and —1), and are a basic tool for the investigation of the 
model at low temperatures. 

In order to study the interface between the two pure phases one needs 
to construct a state describing the coexistence of these phases. Let A be a 
parallelepiped of sides Li,L2,L%, parallel to the axes, and centered at the 
origin of £, and let n = (m, n 2 , n 3 ) be a unit vector in R 3 , such that n 3 ^ 0. 
Introduce the mixed boundary conditions (±,n), for which 



1 if i • n > /p.N 
-1 if % ■ n< U 



These boundary conditions force the system to produce a defect going trans- 
versally through the box A, a big Peierls contour that can be interpreted as 
the microscopic interface. The other defects that appear above and below 
the interface can be described by closed contours inside the pure phases. 

The free energy, per unit area, due to the presence of the interface, is the 
surface tension. It is defined by 

r(n) = hm hm - In (6) 

V ' Li,L 2 ^oo L 3 ^oo /3LiL 2 Z(+)(A) 

In this expression the volume contributions proportional to the free energy 
of the coexisting phases, as well as the boundary effects, cancel, and only the 
contributions to the free energy due to the interface are left. 

Theorem 1. The thermodynamic limit t(vl), of the interfacial free 
energy per unit area, exists, and is a non negative bounded function o/n. Its 
extension by positive homogeneity, /(x) = |x|r(x/|x|) is a convex function 
on R 3 . 

A proof of these statements has been given in ref. [3] using correlation 
inequalities (this being the reason for their validity for all /3) and, in fact, 
the Theorem holds for a large class of lattice systems. The convexity of / is 
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equivalent to the fact that the surface tension r satisfies a thermodynamic 
stability condition called the pyramidal inequality (see [3], jl]). 

Moreover, for the Ising model we know, from Bricmont et al. [5] , Lebowitz 
and Pfister [6J, and the convexity condition, that r(n) is strictly positive for 
T < T c and that it vanishes if T > T c . 

The shape of an equilibrium crystal is obtained, according to the Gibbs 
thermodynamic principle, by minimizing the total surface free energy asso- 
ciated to the crystal-medium interface. The solution to this problem, known 
under the name of Wulff construction, is the following set 

W = {x G R 3 : x • n < r(n) for every n G S 2 } (7) 

Notice that the problem is scale invariant, so that if we solve it for a given 
volume of the crystal phase, we get the solution for other volumes by an 
appropriate scaling. The set W, which will be called the Wulff shape, gives 
the optimal shape for the crystal. 

The Wulff construction can also be viewed as a geometric version of the 
Legendre transformation. Consider the function /(x) defined in Theorem 1. 
From definition ((7|), we get 

W = {x G R d : /*(x) < 0} (8) 

where /* is the Legendre transform 

r(x) = sup(x-y-/(y)) (9) 
y 

Actually /*(x) = 0, if x G W, and /*(x) = oo, otherwise. 

We next introduce a function ip on R 2 such that the graph of x 3 = 
ip(xi,x 2 ), for £3 > 0, coincides with the boundary dW of the crystal shape. 
Since W is a convex body, symmetric with respect to the origin, tp is a 
concave function, and 

W = {x G R 3 : -(f{-xx, -x 2 ) <x 3 < (p(x 1 ,x 2 )} (10) 

This means that — ip is the Legendre transform of the projected surface 
tension r p = (l/n 3 )r, considered as a function on R 2 of the slopes Ui = 
rii/n 3 ,U2 = n2/n 3 . In other words, 



r p (u 1 ,u 2 ) = f(u 1 ,u 2 , 1) 



(11) 
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Indeed, from equations (J5]) and ([9]), we see that 

- <p(u l , u 2 ) = sup {x x ui + x 2 u 2 - r p (u l7 u 2 )) (12) 

U\,U2 

Formula (fT2"j) is known as the Andreev construction [7]). The interest 
of this approach comes from the fact that <p, and hence, the crystal shape 
itself, may be regarded as the free energy associated to a certain statistical 
mechanical Gibbs ensemble. We shall consider this ensemble in sections 2 
and 3. 



2 The two dimensional model 

We next consider the Ising model on a square lattice, with two interaction 
parameters, J\ in the horizontal direction and J 2 in the vertical direction. 
There is in this model an exact expression for the surface tension r(n). How- 
ever, we shall first discuss an approximate expression, which may represent 
this quantity in the low temperature region. It is obtained by restricting the 
sums in formula ()6]) to the ground configurations. 

Let Q°{A) be the set of ground configurations in A under the boundary 
conditions a, and let Af°(A) and E$(A) = H A (a\(±,n)) - H A (+\(+)) be 
the number and the relative energy of such configurations. The set G + {A) 
contains only the configuration (+) and, therefore, 

r(n)= lim lim In V e ^H A [a\(±, a ))+0H A (+](+)) 

= lim lim - — (i^ ±,n) (A) - lnA/" (±,n) (A)) (13) 

Notice that, in dimension two, if n = (— sin 9, cos 9) and t(9) = r(n), 
then v = — tan# and the projected surface tension is r p (v) = r{6)/ cos 9. The 
configurations in (?( ±,n )(A) contain only one Peierls contour, the microscopic 
interface, a polygonal line joining two fixed points in the boundary of A. This 
line is cut only once by all straight lines parallel to the diagonal i\ — i 2 = 0. 
It can then be described by N integers 

0(O),0(l),...,0(iV) (14) 



such that 

n (i) = <p(i) - 4>{i - 1) = ±1, i = l,...,N 



(15) 
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specifying the heights, in units 1 / \/2, over the line %\ + %2 = 0, of the extrem- 
ities of the N consecutive unit segments which compose the polygonal line 
(see Figure 1). Writing u = tan(7r/2 + 9), the boundary conditions are 

(f)(0) = 0, (f)(0) - <f)(N) = uN (16) 



N 




Figure 1. Description of a ground microscopic interface in 2 dimensions. 



Under these boundary conditions we have that n(i) = 1, ^-^N times, and 
n(i) = —1, ^-r^N times. All these configurations have the same energy 

E U (N) = 2J^N + 2 J^N (17) 

and their number is 

wo = (J! N 

2 



Therefore, from equation (TT3T) . we get 

1 N\ 

r p (u) = J\(x+u) + J 2 (i-u) — lim — — In 



N^oofiN (±±±N)\(^N)\ 

. . . . 1 /l+M , 1+u X—u 1 — u\ . . 

J,(i + u) + J 2 (i-u) + - (— In — + — In — ) (19) 
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In order to obtain the equilibrium crystal shape we use the Andreev 
construction (|T2l) . In dimension two, we have 

-</?(£) =svl-p(£u-t p (u)) (20) 



From uq, solution of 



that is 



i = %- = Ji-j*+^— (2i) 

du 2p i-u 



1+U0 = e 2K+2/3.h-2fSJ 2 ^22) 
1— UQ 

e 2/3£+2/3Ji-2/3J 2 _ l 
M ° = e 2/3C+2/3J 1 -2/3J 2 + I ( 23 ) 

we find 

= -£u + Ji + J 2 + u {Ji - J 2 ) + — r ( In h 2 In — - + u In 

Zp * 1 — Mo Z 1 — Uq 

= -£u + J x + J 2 + u {Ji - J 2 ) + (1 + Uo)(£ + Ji - J 2 ) + -g In ^ 
= £ + 2 J x - i In (1 + e 2#+Wi-W 2 ) (24) 
Next, we define the partition function 

e n (£) = r e -W)+w)-*(°)) (25) 

where the sum runs over all configurations with 0(0) = 0. We can interpret 
(j25p as a "grand canonical" partition function with respect to the interface 
boundaries, and the restriction ffl6l) in ffl3l) as a "canonical" constraint. Since 

JV 

0(iV)- 0(0) = $>(z) (26) 
#(0) = 2J x 5(n(z) - 1) + 2 J 2 5(n(i) + 1) (27) 



we have 



(^ = ( e -/3(2J 1 - 5)+e -/3(2J 2 +0)^ ( 2 ^ 



8 



S. Miracle-Sole 



and the corresponding free energy is 

¥>(£)= lim --^ln3jr(£) (29) 
= _ 2_ i n ( e -/J(2*-0 + e -/3(2J 2 +5)) (30) 

= £ + 2 Jx - i In (1 + e 2 ^+ 2 ^~ 2 ^) (31) 

Expressions (T5TT) and (1251) coincide. This means that tp, the free energy 
defined by (1291) . is the Legendre transform of t p , in agreement with the equiv- 
alence of the corresponding "canonical" and "grand canonical" ensembles. 
This fact illustrates the remark at the end of the introduction. 

Putting r] = r}(£) = <p(£), equation (I3T1) can be written 

for 7] > 0. Then, by using the initial coordinates x = 77 + £, y = rj — £, and 
the symmetries of the figure, we obtain the following expression for the Wulff 
shape 

exp(P\x\ - ZK-l) + exp(/3|y| - 2K 2 ) < 1 (33) 

Here Ki = j3Ji, for i = 1,2. This is the approximate solution. We next 
describe the exact solution. 
It has been shown that 

y = y(x) = sup(xv - t p (v)) (34) 

V 

is, up to a factor, the Onsager function for an imaginary argument, 
namely, y(x) = j3~ lA f(—ij3x). This means that 

cosh fly = cosh 2Ki cosh 2K2 — sinh 2K\ sinh IK^ cosh 0x (35) 

where 

K* = -(1/2) IntanhiQ (36) 

is the dual interaction constant (i = 1,2). This is a result due to Abraham 
(see [8]). Taking into account that relation (1361) implies cosh2i^* tanh2Kj = 
1 and sinh 2K* sinh 2Kj = 1, we get 

cosh ^ + cosh ^ - 1 (37) 
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as the exact expression for the Wulff shape. This simple expression has 
not appear in the literature, as far as we know. Several authors (see, for 
instance, [9]) have used the more complicated formula for the surface tension 
r(n), which includes implicit functions, instead of the function y(x). 

The critical inverse temperature (3 C satisfies K\ = K\ or K 2 = K* or, 
equivalently, 

sinh 2K X sinh 2K 2 = 1 (38) 

In the symmetric case, K\ = K 2 = K, we have sinh2iC c = 1 and K c = 
(1/2) \n(V2 + 1) = 0.440687. Now, notice that the left hand side of (ETJ) is 

tanh2if 2 tanl^i^ sinh 27^ + sinh 2^ , . 

> £ _i L = £_! L (39) 

~ cosh2Ki cosh 2K 2 cosh 2K X cosh 2K 2 v ; 

The value of expression (1391) equals 1 for (3 = (3 C , because 

(sinh 2K 2 + sinh 2K x f = 2 + sinh 2 2K X + sinh 2 2K X 

= (1 + sinh 2 21Ti)(l + sinh 2 2K X ) = cosh 2 2K X cosh 2 2K 2 (40) 

This implies that the Wulff shape reduces to the empty set for (3 < f3 c . It 
has a finite positive volume for (3 > (3 C . 

A comparison of expressions ( 1321) and ( 1371) shows that the error made 
when considering the approximate solution (|32|) decreases exponentially in 
f3 when (3 — > oo. However, the shape given by (|32|) presents four corners 
(when the boundary cuts the axes), while that corresponding to ( 1371) has a 
smooth boundary. The same bound on the error can be proven by means of 
a low temperature cluster expansion for the function y, independently of the 
exact solution. This error is less than 1 % for temperatures less than (1/2)T C . 
In the symmetric case the shape defined by (132]) becomes the empty set at 
K' c = (1/2) In 2 = 0.346574. We have K' c < K c . 

3 The three dimensional model: 
Approximate solution 

In this section we consider the Ising model on a cubic lattice and study the 
problem of finding in this case an approximate solution for the equilibrium 
crystal shape, analogous to that of section 2. The system so obtained is very 
interesting in itself, though its validity as a low temperature approximation 
is not known. 
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The ground configurations in C/( ±,n )(A) contain only one Peierls contour, 
the microscopic interface. This contour is a surface made of unit squares 
such that its boundary is a fixed line determined by the boundary conditions 
(±, n) on the faces of the box A. This contour has also the property of 
being cut only once by all straight lines orthogonal to the diagonal plane 
h + ^2 + h = 0. It can therefore be described by specifying the distance, 
or height, of each vertex of the unit squares, which form the contour, to the 
plane %\ + %2 + ^3 = 0. Now, the projection of the cubic lattice Z 3 onto the 
plane i\ + 12 + 23 = is a triangular plane lattice, which will be denoted by 
T. This means that, to each point t E T, we associate an integer variable 
(f>(t) which, multiplied by the factor 1/VS, gives the height of the microscopic 
interface at this point. The height differences 

n(M') = 0(f) - m (41) 

between nearest neighbors t and t' obtained in this way are subject to the 
following restriction: Going clockwise along the edges of each elementary tri- 
angle pointing to the right the values of n(t, t') are {1,1, —2}. An equivalent 
condition is obtained by using the triangles pointing to the left and to follow 
a counterclockwise path (see Figure 2). Conversely, from each configuration 
4>{t), t G T, subject to these restrictions, one can obtain a configuration of 
the microscopic interface, and, hence, a ground configuration in ^^'"^(A), 
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Figure 2. (a) Projection of a cube and heights of the vertices, (b) TISOS 
height differences as obtained from a ground configuration of the triangular 
Ising antifcrromagnet. The corresponding diamond configuration is obtained 
by erasing the edges connecting parallel spins. 



Figure 3. Illustration of the mapping of a ground microscopic interface in 
3 dimensions onto a diamond configuration. A diamond configuration is 
obtained by erasing one edge in every elementary triangle of the triangular 
lattice. 
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provided that the boundary conditions are satisfied. The model with such 
a set of configurations was introduced in ref. [TJj and considered also in ref. 
[T2] . It is called the tisos (triangular Ising solid-on-solid) model. 

There is a one-to-one correspondence between the tisos configurations 
and the allowed configurations of the Ising antiferromagnet at zero temper- 
ature. One can obtain it by requiring that odd and even height variables 
correspond to spins of different sign. It is easily checked that this Ising con- 
figuration obeys the T = antiferromagnetic constraint that no elementary 
triangle should contain three spins with equal signs. 

If one erases all lattice edges between parallel spins, one obtains a plane 
filled with rhombi or diamonds, as shown in Figure 3. Such a diamond 
covering is just the orthogonal projection on the diagonal plane of the edges 
of the microscopic interface initially considered on the cubic lattice. It gives 
also a view in perspective of this interface in the three dimensional space. 

Let ei, e 2 and e3 be three vectors on the plane at angles 27r/3 and modulus 
equal to y/3/2. Since the length of the edges of T is a/3/2, the sites of the 
lattice T are of the form + t 2 e 2 , with integer coordinates ti,t 2 . Let 
Q be the set of lattice sites inside the parallelogram defined by \t\\ < L\ 
and \t 2 \ < L 2 . Let \Q\ be the area of this parallelogram and let dQ be the 
boundary of Q (i.e., the set of sites in Q having some neighbor outside Q). 

The boundary conditions corresponding to the slope u = (ui,u 2 ) are 

<f>(t) = <j> u (t) = [mi! + u 2 i 2 ] } % e dQ (42) 

where [•] denotes the integer part. The corresponding projected surface ten- 
sion is given by 

t p (u 1 ,u 2 ) = lim —oprrZ A (Ui,U 2 ) (43) 
Li,L 2 -s-oo p\Q\ 

with 

Z A ( Ul ,u 2 ) = Y,e- m<p) (44) 

In ( 1441) the sum runs over all configurations inside Q satisfying the bound- 
ary conditions ( )42l) . and H(4>) = i?A (<?"</> n )) — -^a(+I(+)) is the relative 
energy of the ground configuration associated to the microscopic interface 
defined by 0. 

We can interpret the conditions ( 1421) as "canonical" constraints and intro- 
duce a conjugate Gibbs ensemble of ( 1441 . which can be viewed as a "grand 
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canonical" ensemble with respect to the interface boundaries. For this pur- 
pose we consider the boundary terms 

SM= E S M = E (45) 

£e£i(Q) tet 2 {Q) 

where £%(Q) and £ 2 (Q) are the sets of all bonds in Q parallel to the vectors ei 
and e 2 , respectively, and oriented according to increasing coordinates. The 
bond variables n(£) are the height differences of formula f!4ip . The grand 
canonical prescription, which is convenient to consider, consists in adding to 
the energy a term of the form 

xxS x ((f>) + x 2 S 2 (<f ) ) (46) 

where x = (xi,x 2 ) G R 2 represent the slope "chemical" potentials. The 
associated partition function and free energy are 

E p Q er {x 1 ,x 2 ) = ]T e-WM+PxiSi M+Px 2 s 2 w (47) 
4> 

The sum in (j47p runs over all configurations in Q with periodic boundary 
conditions (with respect to the bond variables n(£)), 0(0) is taken equal to 
0. 

Theorem 2. The limit, when Q — > oo, of tpq 1 , exists and equals 

(p(xi, x 2 ) = - sup{a;iMi + x 2 u 2 - t p (ui, u 2 )} (49) 

u 

in the interior of the domain of the concave function if. 

The theorem is proved in ref. [10] for many sos models. Notice that the 
partition function (j4"81) does not correspond exactly to the conjugate ensemble 
of the function (I44p . as conditions (1421) represent many more constraints that 
just fixing the values of the two sums in (j450 . The domain of a concave 
function <p(x) is the set {x G R 2 : <f(x) > — oo}. 

It turns out that the free energy function ip(x) of the tisos model can 
be computed exactly. There is no additional difficulty in considering, like in 
section 1, an asymmetric Ising model on the cubic lattice (with three different 
interaction parameters in the vertical and the two horizontal directions). 
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However, it is interesting instead to introduce a new kind of asymmetry, 
and to take the interaction parameters equal to J\, J 2 or J 3 , according to 
whether the coordinates of the highest point of the corresponding bond satisfy 
H + ^2 + ^3 = 0, 1, or 2 (mod 3). 

Then, with Ki = {3Ji for i = 1,2, 3, we write 

n = exp (-Ki + K 2 + K 3 ) 
r 2 = exp (K ± -K 2 + K 3 ) 
r 3 = exp (Kt + K 2 - K 3 ) 

r 3 = n r 2 r 3 = exp (K x + K 2 + K 3 ) (50) 



and 



R = (ri + 4 + 4)r i (51) 



For reasons of symmetry, the grand canonical prescription, which will be 
considered, consists in adding to the energy, instead of the term (146]) . a term 
of the form 

xiS 1 ((f>)+x 2 S 2 ((f>)+x 3 S a ((f>), withS 3 (0)= ]T n(£) (52) 

ie£ 3 (Q) 

where x 3 G R and £ 3 {Q) is the set of all bonds in Q parallel to the vector e 3 . 
Taking, X\ = (3xi, for i = 1, 2, 3, we write 

w x = exp (-Xl + X 2 + X 3 ) 
w 2 = exp (Xi - X 2 + X 3 ) 
w 3 = exp (Xi + X 2 - X 3 ) 

Wo = wiw 2 w 3 = exp (Xi + X 2 + X 3 ) (53) 

Then, we have 

2 _ 

93 = lim — — In 5 = ln(w /r ) + 

il,i2->00 L\L 2 

— ^ / 2 " rfp / 2 " In litojje* + w 3 e 3i P _ w 3 e -3ig + ^ 
247?^ jo jo 

The exact solution (15~4"|) for the free energy of the model is due to Nienhuis 
et al. [12] (where slightly different notations are used). A careful study of 
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this solution, also in ref. [12], shows the existence of four regions in which 
the function (p has simple expressions, as indicated below 

1) if —w\ + wl + wf + Rwq < 0, <p = Inwx + In (w /r ) 

2) if w\ — wl + w| + Rwq < 0, = l nu; 2 + In (w /r ) 

3) if tyf + tyf — u>| + -R Wq < 0, = In w; 3 + In (w /r ) 

4) if w\ + tof + w| — Rwq < 0, = y3 



In regions 1), 2) and 3), <£> is a linear function of the variables x\,X2,X3. 
It is constant in region 4), where it takes the value of fl54|) at the point 
x\ = %2 = x% = 0. These regions correspond to the domains of four ordered 
phases of the tisos model in the phase diagram of this system. 

According to the Andreev construction, the graph of the function 

z = -{l/P)<p{xi,x 2 ) (55) 

gives the boundary of the Wulff shape for the equilibrium crystal (equation 
fflQ]) ). The factor —1//3 comes from definition fl54l) in order to obtain the 
correct free energy. We have, in region 3), for instance, 

z = h + h + Js - 2(xi + x 2 ) (56) 

which shows that the Wulff shape contains a plane facet orthogonal to the 
(0,0,1) direction. The boundary of this facet is determined by the boundary 
of region 3). The other regions 1), 2) and 4) indicate the existence of three 
other facets respectively orthogonal to the (1,0,0), (0,1,0) and (1,1,1) direc- 
tions. The portion outside these regions corresponds to the curved part of 
the crystal surface. 

It is only in recent times that equilibrium crystals have been produced in 
the laboratory, information and references on this subject can be found in 
the review articles pQ and [2]. As in the present model, a typical equilibrium 
crystal at low temperatures has smooth plane facets linked by rounded edges 
and corners. The area of a particular facet decreases as the temperature is 
raised and the facet finally disappears at a temperature characteristic of its 
orientation, a phenomenon known under the name of roughening transition. 
The reader will find in ref. [12] a study of the roughening transition for the 
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(1,1,1) facet in the tisos model, as well as a discussion on several other 
interesting properties of the crystal shape. 

In the symmetric case, in which the interaction parameters of the Ising 
model are equal, J\ = J2 = J 3 = J, we have R = 3 (from equations (150]) and 
( JBTI) ). Then the region of phase 4) is empty and the (1,1,1) facet does not 
exist. This follows from the algebraic identity 

w l + w 2 + w 3 ~ 3W1W2W3 = 

(Wi +W 2 + W 3 )(wf + w\ + wl - WiW 2 - W 2 W 3 - W 3 Wi) = 

+ W 2 + W 3 )((W! - w 2 f + (w 2 - w 3 ) 2 + (w 3 - Wif ) (57) 

showing that 

w\ + w\ + w\ - 3w% > (58) 

Conversely, in the case that J 1; J 2 and J 3 are not equal one has R > 
3 and expression f )5^]l can be negative. Then region 4) is not empty if 
is large enough, which implies the appearance of the (1,1,1) facet at low 
temperatures. 




Figure 4. Equilibrium crystal shape according to the TISOS model in the 
symmetric case. The crystal is shown in a projection parallel to the (1,1,1) 
direction. The three regions 1), 2) and 3) indicate the facets, and the re- 
maining area represents the curved part of the crystal surface. 
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Figure 4 shows a (1,1,1) corner of the crystal for the symmetric case 
J\ = J2 = J3 = J- In this case the definition of region 3), for instance, can 
be written, using again identity f l57|) . 

wf + w\ — wl + 3W1W2W3 = 

+ w 2 - w 3 )(( Wl - w 2 ) 2 + (w 2 + w 3 ) 2 + (-w 3 - Wl f) < (59) 

or, equivalently, 

w\ + W2 — w 3 < (60) 
and, using fl53|) with X3 = 0, we see that region 3) can be defined by 

The comparison of this formula with equation (I3"3"|) shows the following state- 
ment. 

Remark. The shape of the facets of the equilibrium crystal associated 
with the tisos model, in the symmetric case, coincide with the Wulff shape 
obtained from the approximate solution of the two dimensional Ising model. 

This property can be derived, without using the exact solution (j5i|) . with 
the help of the results of section 4 below, in particular Corollary 1. These 
results apply also to the tisos model, and the shape of the facet is obtained 
by computing the corresponding r step (m). Notice also that the model can be 
used to describe the (1,1,1) corner of the crystal, that is 1/8 of the crystal. 
Completing it by symmetry we see that its facets undergo a roughening 
transition at K' c = (1/2) In 2 = 0.346574 (see last paragraph of section 2), 
and we have K' c > K c (3), where K c (3) ~ 0.22 is the value of the critical 
coupling constant for the three dimensional Ising model. 

There are other models for which an analysis similar to that of the present 
section could be developed. For example, the Ising model on a body centered 
cubic lattice with nearest and next nearest neighbor interactions is one of 
these models. This system leads also to an approximate solution that can be 
expressed in terms of exactly solvable sos models (see [13] and the references 
quoted there). 
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4 The three dimensional model: 
Rigorous results 

In the previous section we have described the microscopic interfaces at zero 
temperature for the Ising model on a cubic lattice. Some properties of these 
interfaces were studied by means of the exact solution of an associated sos 
model. At low (positive) temperature we expect the ground interfaces to be 
modified by deformations. Small deformations will appear here and there, 
the large deformations having a very small probability. It appears, however, 
very difficult to develop such an argument into a rigorous proof that could 
justify the validity at low temperatures of the approximate solution described 
in section 3. 

The situation is different for the interface orthogonal to the direction 
n = (0,0,1). In this case there is, at zero temperature, only one micro- 
scopic interface, which coincides with the horizontal plane. The small defor- 
mations of this interface that appear at low temperatures are accessible to 
a mathematical treatment. The interface is rigid at low temperatures and 
its properties can be studied by means of a convergent cluster expansion. 
By pursuing this analysis it has been possible to determine the shape of the 
facets in a rigorous way. It can then be seen [13], from the proof of Theo- 
rem 3 below and Corollary 1, that this shape differs from the Wulff shape 
associated with the two dimensional Ising model, only by a quantity of or- 
der exp(— 12/3 J), and thus, that it is close to the shape predicted by the 
approximate solution of section 3. 

The appearance of a facet in the equilibrium crystal shape is related, 
according to the Wulff construction, to the existence of a discontinuity in 
the derivative of the surface tension with respect to the orientation. More 
precisely, assume that the surface tension satisfies the convexity condition of 
Theorem 1, and let this function r(n) = t(9, (ft) be expressed in terms of the 
spherical co-ordinates < 9 < n, < <p < 2n of n, the vector no being taken 
as the X3 axis. A facet orthogonal to the direction n appears in the Wulff 
shape W if, and only if, the derivative dr(9,<j))/d9 is discontinuous at the 
point 9 = 0, for all 0. Moreover, the one-sided derivatives with respect to 9 
exist, at 9 = + and 9 = - , and determine the shape of the facet. The facet 
J 7 C dW consists of the points x G R 3 belonging to the plane £3 = T(no) 
and such that 

Xi cos (ft + x 2 sin (p < dr{9, <f))/d9\ g =0 + (62) 
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for every between and 2tt. A proof of this fact can be found in ref. [H] 
(Theorem 1). 

The step free energy is expected to play an important role in the facet 
formation. It is defined as the free energy (per unit length) associated with 
the introduction of a step of height 1 on the interface, and can be regarded 
as an order parameter for the roughening transition. Let A be a paral- 
lelepiped of sides Li,L 2 ,L 3 , parallel to the axes, and centered at the origin, 
and introduce the ( step, in ) boundary conditions, associated to the unit vec- 
tors m = (cos 0, sin 0) G R 2 , by 

-/■\_ fl if i > or if i 3 = and i\m\ + i^rriz > 
I —1 otherwise 

Then, the step free energy, for a step orthogonal to m (with m2 7^ 0) on the 
interface orthogonal to n = (0, 0, 1), is 

r stcp (0) = lim lim lim In ^ w l } (64) 

A first result concerning this point, was obtained by Bricmont et al. [15J. 
These authors proved a correlation inequality which establish r step (0) as a 
lower bound to the one-sided derivative dr(9, 0)/d9 at 9 = + (the inequality 
extends also to 7^ 0, see [H]). Thus when r step > a facet is expected. 

Using the perturbation theory of the horizontal interface, it is possible, as 
mentioned above, to study also the microscopic interfaces associated with the 
(step, m) boundary conditions. The step structure at low temperatures can 
then be analyzed with the help of a new cluster expansion. As a consequence 
of this analysis we have the following theorem. 

Theorem 3. If the temperature is low enough, i.e., if (3 > /3q, where 
0o > is a given constant, then the step free energy r step (0), defined by 
limit (64\ ), exists, is strictly positive, and extends by positive homogeneity to 
a strictly convex function. Moreover, r stcp (0) can be expressed as an analytic 
function of (5, for which a convergent series expansion can be found. 

Using the above results on the step structure, similar methods allow us 
to evaluate the increment in surface tension of an interface titled by a very 
small angle 9 with respect to the rigid horizontal interface. This increment 
can be expressed in terms of the step free energy and one obtains the following 
relation. 
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Theorem 4. For f3 > (3q, we have 



dT(9,(f>)/d6\e =0+ = r 



(65) 



This relation, together with equation (|62|) . implies that one obtains the 
shape of the facet by means of the two-dimensional Wulff construction applied 
to the step free energy. 

Corollary 1. The equilibrium crystal presents, if (3 > 0o, a facet 
orthogonal to the (0,0,1) direction, whose shape is given by 



The reader will find a detailed discussion on these points, as well as the 
proofs of Theorems 3 and 4, in ref. [H]. A brief explanation of these results 
can be found in ref. [16]. 

From the properties of r step stated in Theorem 3, it follows that the 
Wulff equilibrium crystal presents well defined boundary lines, smooth and 
without straight segments, between a rounded part of the crystal surface and 
the facets parallel to the three main lattice planes. 

It is expected, but not proved, that at a higher temperature, but before 
reaching the critical temperature, the facets associated with the Ising model 
undergo a roughening transition. The conjectured value of /3r, the roughen- 
ing inverse temperature, is given by Kr ~ 0.38 [17J. It is then natural to 
believe that the equality (j65|) is true for any (3 larger than /3r, allowing us to 
determine the facet shape from Corollary 1, and that for < (3r, both sides 
in the equality vanish, and thus, the disappearance of the facet is involved. 
However, the condition that the temperature is low enough is needed in the 
proofs of the above results. 
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